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1 Introduction 



It is of interest to search for 5d gravity duals [lj of time dependent phenomena in 4d gauge 
theories [21 [31 01 [6J . These could, for example, serve as prototypes of the dynamics of heavy 
ion collisions. For quasi-static phenomena in strongly coupled matter gauge-gravity duality 
picture has already produced lots of interesting results, for viscosity [7J [U [9], jet energy loss 
PHI QH CE21 UJ] and for photon production [14]. 

One can model collisions of large nuclei at very high energies by taking the transverse size 
and collision energy to be effectively infinite so that the dynamics will be invariant under 



longitudinal (xi) Lorentz transformations. Natural variables then are r = \jt 2 — x 2 , r] = 
l/21og(xi/£) and there will be no dependence on X2,x%. The 5d gravity dual of this 4d 
physics then has the metric [2] 



ds 



r 2 - 
2 = — 

z 2 



-a(r, z)dr 2 + r 2 6(r, z)drj 2 + c(r, z){dx\ + dx\) + dz z 



The unknown functions a, b and c would be determined as solutions of 5d AdS gravity equa- 
tions and 4d physics on the boundary at z = could be computed. Several interesting and 
suggestive results have been obtained by studying large-r behavior, without using an exact 
solution. In particular, perfect fluid behavior for the energy density e(r) ~ 1/r 4 / 3 is singled 
out by absence of curvature singularities [2], also evidence for the viscosity /entropy density 
ratio r]/s = h/An has been found by studying corrections to the leading large-r behavior [6j. 

To clarify the issue we will here restore the 3d spherical symmetry broken by the metric 
(TjQ). For this case we found an exact solution, which, in fact, proved to be a known solution 
[151 [TBI E] i n different coordinates. The boundary theory is the isotropic and homogeneous 
cosmological FRW metric with an arbitrary scale factor r{t), it is thus rather big than little 
bang. In big bang energy density decreases because space expands. Here we have no 4d 
gravity to determine r(t), but we shall fix it so that the comoving energy density decreases as 
that in the center of a spherical little bang. Little bang takes place in Minkowski space and 
the energy density in the center, in the rest frame, decreases because matter flows outwards. 
To permit one to go away from the rest frame or the comoving frame would require a more 
complicated ansatz for the 5d metric than what we use in this paper. 



2 Similarity expansion in relativistic hydrodynamics 

Experimental measurements of QCD matter in relativistic heavy ion collisions have shown 
that the matter flows like a nearly ideal fluid. As discussed above, for large nuclei and ener- 
gies the expansion can be approximated by a 1+1 dimensional longitudinal boost invariant 
similarity flow of 3-dimensionally thermalised matter, 

v i (t,x\x 2 ,x 3 ) = ?j5 il . i = 1,2,3. (2) 

Let us now generalise this to radial flow in d — 1 spatial dimensions, keeping d arbitrary for 
the moment. 
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A radial similarity flow is defined by the velocity field 
v(i,x) = y, n^ = (7, 7 v) = ^=, M = 0,l,..,d-1, -x 2 = -rj^x "x v = t 2 - x 2 = r 2 



(3) 

The energy-momentum tensor is 

2}^ = (e + p)u^u u +p<q ilv = {e + p) + p rj^ (4) 

and the conservation equation d^T^ v = 0, split in components parallel and perpendicular to 
u^, can be written in the form 

x"<V + (d-l)(e + p) = 0, (5) 
^ V -—)d»P = 0. (6) 

The second is solved by p(t, x) = p{r) and, if the conformally invariant equation of state 
p = e/(d — 1) is assumed, the first one becomes 

re'(r) + de(r) = =* e(r) = g = ^ _ £ ° y/2 . (7) 

The constant eo here is given by the initial conditions of the flow. Formally, it could also be 
fixed by including a source term on the RHS of ([5]). 

At any fixed time t the flow thus consists of a spherical shell with infinite energy density 
moving radially with light velocity. In the interior the flow pattern is of the similarity type 
([3]) with energy density only depending on the proper time. In the local rest frame x = 

e(r) = e(t) = |. (8) 

Our goal thus is to find the gravity dual of radially expanding thermalised matter with the 
energy-momentum tensor Q and comoving energy density dH). 

One may try to express T^ u in terms of radial rapidity coordinates (T,rj r ,8,(p), t = 
Tcosh.r] r , r = rsinhr/ r ; 9,(j) are the usual spherical angles. Then x M — > (r, 0,0,0) but 
Viw —* 9nu with complicated non-diagonal terms. 

In the above matter is flowing in Minkowski space; compare it now with space expanding 
as m the flat FEW metric diag(-l, r 2 (i), ..,r 2 (t)). Then 

(d - l)r'(i) , » , x 

£(f)+ r(t) ( £ + P) = - ( 9 ) 

One observes that if r(t) ~ t this coincides with (0) evaluated in the local rest frame v = 
0, x = 0. For this behavior of the radius the energy density decreases in the same way for 
radial flow in Minkowski space and for Hubble flow in the FRW metric. We shall later in 
Section [5] show how gauge-gravity duality selects just this power. This is in analogy to the 
selection [2] of the adiabatic power e ~ 

]y T 4/3 for the flow ^ 
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For completeness, the standard dissipative tensor for the flow ([3]) is 



AT - = - c 7^("--^)- (10) 

The shear part of the dissipative tensor vanishes and thus studying this flow does not give 
a handle on the shear viscosity. Eq.([7j) becomes re'(r) + de(r) — Q{d — l) 2 /r = 0, but for a 
massless conformal fluid further £ = 0. 



3 The exact solution 

In this article we shall search for gravity duals of spherically expanding systems; the perfect 
fluid solution was presented in detail in the previous section. The gravity solution thus has 
to be of the type 



ds 



r 2 

2 L - 



-a(t, z)dt 2 + b(t, z)(dx\ + dx\ + dx\) + dz 2 



C 



2 



= guNdx M dx N = gu U (x, z)dx^dx v + dz 2 , (11) 
z z L J 

where the metric is determined from the 5d AdS gravity equations 

T^MN ~ 2^9MN ~ -j^9MN = 0, x M = (t, x 1 , x 2 , x 3 ,z). (12) 

Here the AdS radius C is related to the 5d Newton's constant by 

£ 3 2N 2 



G fi " vr 



(13) 



The functions a(t, z) and b(t, z) in (jlip can now be solved from (|12p as follows |18j . Using 
the tz-component of (|12() . one can solve a(t,z) in terms of an arbitrary function Fi(t): 

b(t,z) 

Using this solution for a(t, z), one can solve the function b(z, t) from the it-equation in terms 
of two more arbitrary functions: 



Kt,z)=F 2 (t)z 2 + ^l-^ ry (15) 



Finally, using the zz-equation one gets a differential equation for iq(i), F%(t) and Fs(t), 
from which one of the functions can be eliminated. We specify the remaining two arbitrary 
functions r(t) and h(t) so that the 4d boundary metric is 

g$ = di ag (-h 2 (t),r 2 (t) ir 2 (t),r 2 (t)) (16) 
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and the constant of integration \/2zq so that it is the position of the horizon in z when 
r = h = 1. One can trivially set h = 1 by a choice of time coordinate and the boundary 
metric then is a flat cosmological FRW metric. 

The solution is, abbreviating r = r(t), h = h{t),r' = r'(t),r" = r"(t), h! = h'(t), 



2J1 



a(t, z) 



h r 
b(t,z) 



1 + 



/ h'r' 



\2h 3 r 2h 2 r 
+ 



z 2 + 



„/4 



h'r 



8/i 4 r 3 ~ 16/i 4 r 4 ~ 8h 5 r 3 ~ 4r 4 4 



h=l 



b(t,z) 



Ar 



b(t,z) 



„'2 



4/ l 2 r 2 



Z* + 



l 



4r 4 ; 



4r 4 ; 



o 



For r(t) = h(t) = 1 the solution reduces to the metric [2] 



ds 2 = 



(l-z 4 /(4z 4 )) 2 



1 + z 4 /(44) 



dt 2 + 1 + — -j\ (dx{ + dx 2 2 + dxg) + dz' 



4z 4 



(17) 
(18) 



(19) 



which, choosing a new variable z 2 = z 2 /(l + z 4 /4zq), can be brought to the standard AdS 
black hole form 



ds* 



z 2 



z 

-(1 -£)dt 2 + rfx? + dx\ + cfa| + 



1 



1 - z 4 /z 4 



(20) 



The metrics have a horizon at z = V2zq,z = zq, with Hawking temperature Tjj = 1/(itzq) 
and with entropy density A/(AG 5 V 3 ) = C 3 /{AG 5 z^) = tt 2 N 2 T^/2. 



4 Energy-momentum tensor 

To study the dynamics of the boundary theory, we have to find its covariantly conserved 
energy-momentum tensor. We present two methods for doing this. 

For the first [19], expand the g^ u in (fTTj) . (fT71) and (fl~8j) near the boundary: 

9At,z)=g$(t) + g$(t)z 2 + g$(t> 4 + ---, (21) 

where g^ is the FRW metric in (|16|) and the rest are easy to work out from (I17p and (I18p . 
Then the energy-momentum tensor can be evaluated using 



» v ~ 4vrG 5 



(22) 



" |<W 2 >) 2 - Tr (gM) 2 } - \(g^ g~^ 2 \ v + \(Tr g®) • g® 
with the result 

Tt = g^T au = diag(-e(t),T 1 1 (t),T 1 1 (i),r i 1 (t)), (23) 
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where, using (fT3|) . 



m1 , N 1 , . iV r 2 r n (h'r'-hr") , s 



i 37V 2 r l2 (h'r' - hr") , s 



The trace of is 



which is just the standard trace anomaly [20J. Further, the tensor T^ v is covariantly conserved, 
WT^„ = 0, leading f or v = to 

e '(t) + —[ e (t)+r 1 1 (i)] = 0. (27) 

For the second, the same result can be obtained without the expansion of g^ u (x,z) by 
noting that, (a) due to the imposed symmetries must have the form (|23j) . (b) must 
be conserved as in (p?]) . (c) the trace anomaly is C 3 /(64:TrG 5 )(R flu R> MU - R 2 /3), which leads 
to (|26p . Eliminating using the anomaly equation leads to an equation for e(t) which can 
be solved to again give the result in (|2^f) and (f25|) . 

The boundary energy-momentum tensor can be naturally interpreted to describe mass- 
less fluid in a curved background metric. The energy density has two components: e(t) = 
eo(i) + Ae(i). The temperature dependent first part eo(t) ~ r(t)~ 4 is the standard behav- 
ior of homogeneous radiation in expanding spacetime and the temperature independent part 
Ae(t) ~ r'(i) 4 /r(£) 4 describes quantum corrections to the matter due to the curved back- 
ground. 

Also the Tij can be suggestively composed to two parts: Tjj = p(t)S{j — ATij, where 
p(t) = 1/3 e(t) is the pressure of the fluid. One may try to interpret AT™ in terms of 
bulk viscosity: comparing (|25|) with the dissipative part (Eq. (|10p for x = 0) of the energy- 
momentum tensor, AT^ = — CfcV • v = — 3£r' /rSu we can identify a time dependent bulk 
viscosity (take h = 1) 

iV? rV , . 

cw = ^— • m 

Bulk viscosity vanishes for a conformal massless system, but here it is precisely the anomaly 
which gives rise to it. A positive £ implies that entropy increases as can be explicitly verified 
from our solution. 

5 Entropy, gravity dual 

The coefficient a(t, z) in (|17p can be written in the form (for h = 1) 

„2 



r 

a(t, z) = — r 

V ' b(t, z) 



1 




2 

(29) 
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where 



Z H± 



Ar 2 



ll\2 



ZH+ = Z H {t) < Z H - 



(30) 



rr" ± y^/z 4 + (r' 2 

We have now a bulk theory with a dynamical horizon at zji(t). The entropy density s(t) then 
is given by the area of the horizon: 



s(t)r 3 (t) = £ 
^3 



Area 



1 



AG 5 V 3 4G 5 



d 3 x 



C 3 b 3/2 N 2 b 3/2 ( t ,Z H ) 



2-Kzjj 



where 7 is the determinant of the metric of the (a; 1 ,^ 2 ,^; 3 



b(t, z H ) _ 4/4 + K' - r ' 2 + \j^/4 + i r ' 2 - rr ") 2 } 



4G 5 z 3 H ' 

subspace in (fTT|) and 

2 



A(rr" + a/4/4 + (r 



12 



"\2\ 



(31) 



(32) 



For arbitrary r(t) we thus have the entropy, but the issue of defining a temperature for a 
dynamical horizon is a very complicated one [21J. It should be defined so as to satisfy the 
thermodynamic relations e+p = Ts, s(T) = dp/dT, for thermal energy density and pressure. 

Further, it is of interest to compute the curvature invariants for the solution (|18p in order 
to study whether the bulk solution is singular or not [2j. Because (|lip is a solution to 
the equations of motion, we automatically have 1Z = —20/C 2 , 1Z mn 1Zmn = 80/£ 4 . The 
expression for 1Z M np ®1Zmnpq simplifies to the form 



n MNp Qn MNPQ = -L \ 40 + 72 



zgb(t, z) 



(33) 



where b(t, z) is given by (|18p . The 40 here is the maximally symmetric part 21Z 2 /(d 2 — d) of 



'*-M NPQ ■ 



To calibrate the above expressions note that for r = 1 we have z 2 u = 2zq and b{zH)/z 2 j = 
1/z 2 , which, using (f3Tj) . gives s = N 2 /(2-ttzq), the standard value for a static horizon. Fur- 
thermore, the metric then is (1201) with the Hawking temperature T = 1/{ttzq) and one obtains 
the standard result s = tt 2 N 2 T 3 /2,p = ir 2 N 2 T 4 /8 for the pressure of strongly coupled ther- 
malised super symmetric Yang-Mills matter. Also, for r = 1 the curvature invariant (1331) has 
the value T1 2 mnpq = H2/£ 4 at the horizon. 

Our goal is to find the gravity dual of the flow in ([4]) with a powerlike e(t) in ([8]) with 
d = 4. For a general powerlike behavior r = (t/to) n we have r' = nr/t, r" = n{n — l)r/t 2 so 
that 




z y 



4f 2 




3n 4 + 4ra 3 



n 2 -n+ ^/4t 4 /(z 4 r 4 ^ 
4i 4 /(z 4 r 4 



+ n 2 



Aira—^- 



+ 



-n + ^/4£ 4 /(z 4 r 4 ^ 



n 2 - n + yj^/izftr 4 ) + n 2 



3/2 



(34) 
(35) 

(36) 
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where a = N 2 /(8ir 2 ). All the quantities depend essentially on the combination 4t 4 / (z 4 r 4 (f)) . 
The required behavior e = e^/t 4 is thus obtained when n = 1, r(t) = t/tQ. Then 

2t 3a { , *\ 3n 2 N 2 f y/2 



ZH = (4^ + l)V4. e =4FU + 1 J = ^l^)J =3P ' (37) 




. n \ 3/2 

2vra / /4tS \ ' n*N? I I a/2 \ 1 



t 3 1 V 4 / 2 I \TTZ H (t) I 2irH 2 



(38) 



The gravity dual we searched for thus is given by the metric (|17p . (|18p with r(t) = t/tQ. 
This result has the following properties: 

• For the boundary flow (jHJ) the scale factor eo is given by initial conditions. Here this is 
related to the parameter to of the dual metric. Writing 1/zq = ttTjj it appears in the 
combination AttHqT^ + 1, where the factor +1 represents the effect of the curvature of 
the boundary metric. We want this to be negligible, since the little bang takes place in 
a flat space. Thus we have to demand 

7TT H t = - > 1. (39) 

Then also the scale factor in (jSJ has to satisfy eo > 3a = 3N 2 /(8ir 2 ). 

• Since we are looking for the gravity dual of a thermalised system, we also have to define 
the temperature. This is defined by the two conditions T = (e + p)/s = Ap/s and 
s = dp/dT. These can be integrated to give p = cT 4 , s = 4cT s , where c = (s/4) 4 /p 3 
is a constant, given by Eqs. (|38|) . For the temperature one obtains from T = 4p/s and 
Eqs. d3HJ that 

4 4j.4t->4 1 1 

v / 2vrrt = 5_JLZ ^V^Ttfio ii irT H t >l. (40) 

(y/4**4T& + l-l) 

Inserting this to s = 4cT 3 , p = cT 4 , one sees that the equation of state s(t) = 
7} ir 2 N 2 T 3 (t) and p(t) = ^iT 2 N 2 T 4 (t), appropriate for the matter we are consider- 
ing, is obtained if irTnto = to/zo > 1. Thus matter with correct equation of state at a 
temperature satisfying Tt = Tnto is defined in the boundary theory when (|39|) holds. 
The behaviour Tt ~ constant is the analogue of the Tt 1 / 3 = constant behaviour for a 
longitudinally expanding system. 

• Eq. (|39|) is also a natural condition for the initial temperature of a matter system 
thermalised at to a t the initial temperature Th, just from the uncertainty principle. 
For to<zo = 1/ttTh one is in a pre-thermalisation region in which T, e = 3p and s are 
formally defined by the above equations but are not related by the correct equation of 
state. In this pre-thermalisation region curvature effects are large. 
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• The adiabatic power e ~ 1/t 4 / 3 for the flow in ([2]) was singled out in [2] by demanding 
regularity of the invariant TZ?mnpq a ^ the horizon. The same argument works out 
here, too. Evaluating (|33p at the horizon for r(t) ~ t n shows that it is constant for 
n = 1, grows without bounds ~ t 8 ( n_1 ) for t > when n > 1 and grows without 
bounds ~ l/t 8 ( 1_n ) for t <C zq when n < 1. Thus, if one demands that T^mnpq & t the 
horizon be regular when t — > oo or t — ► 0, one again finds that n = 1, corresponding to 
e(i) = eo/t 4 , is the only option. 

6 Conclusions 

The goal of this work was to develop the framework for studying time dependent systems 
in the gauge/gravity duality picture in [2J [31 [H [6] by analysing an exact solution of the 
5d AdS gravity equation. The price we had to pay was that more symmetry had to be 
assumed: we had as the goal 3d radial expansion. It turns out that the boundary metric 
is the standard cosmological FRW metric with an unknown time dependent radial function 
r(t). In cosmology one, of course, can determine r(t) from Einstein's equations when T^ v 
is given. Since we want to have adiabatically expanding matter in the boundary theory, we 
determine r(t) = t/to so that the correct T^ u for the expanding system is obtained. The 
uncertainty in to corresponds to arbitrariness of the initial conditions of the radial flow, for 
toTo>l the system is thermalised and curvature effects of the boundary metric are negligible. 
Conversely, by demanding regularity of T1 2 mnpq & t the horizon, one can derive the adiabatic 
exponent in e ~ 1/i 4 . 

It would be most interesting to also obtain information on the transport coefficients. How- 
ever, with spherical symmetry, v = r/t, shear viscosity does not contribute, Eq. (jlOh . and 
for the perfect fluid case r(t) ~ t the anomaly, interpreted as bulk viscosity, vanishes. 

There are clearly many issues requiring further study. For the first, an exact solution with 
the symmetries ([TJ would be very useful. The boundary metric would then be of the form 
diag(— 1, g 2 (t) , r 2 (t) , r 2 (t)) [5\. For the second, consequences of the time dependent horizon 
and its effects on the determination of entropy and temperature should be understood better. 
The passage from known energy density and pressure, with non-thermal components, and 
entropy to temperature is the key issue here. For the third, exact solutions with scalar and 
form fields could possibly give more information. 

Acknowledgements. We thank Esko Keski-Vakkuri for discussions and K. Skenderis, J. 
Louko and R. Janik for advice. This research has been supported by Academy of Finland, 
contract number 109720. 

References 

[1] J. M. Maldacena, "The large N limit of superconformal field theories and supergrav- 
ity ," Adv. Theor. Mat h. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 (1999)] 
arXiv:hep-th/9711200] . 



8 



R. A. Janik and R. Peschanski, "Asymptotic perfect fluid dynamics as a consequence of 
AdS/CFT," Phys. Rev. D 73, 045013 (2006) EarXiv:hep-th/0512162] . 

R. A. Janik and R. Peschanski, "Gauge / gravity duality and thermalization of a boost- 
invariant perfect fluid," Phys. Rev. D 74, 046007 (2006) | arXiv:hep-th/0606149] . 

S. Nakamura and S. J. Sin, "A holographic dual of hydrodynamics," JHEP 0609, 020 
(2006) [arXiv:hep-th/0607123| . 

S. J. Sin, S. Nakamura and S. P. Kim, "Elliptic flow, Kasner universe and holographic 



dual of RHIC fireball," |arXiv:hep-th/0610113 



R. A. Janik, "Viscosity from nonsing ularity," |arXiv:hep-th/0610144| 

G. Policastro, D. T. Son and A. O. Starinets, "The shear viscosity of strongly cou- 
pled N = 4 supersymmetric Yang-Mills plasma," Phys. Rev. Lett. 87, 081601 (2001) 
[arXiv:hep-th/0104066l . 

P. Kovtun, D. T. Son and A. O. Starinets, "Viscosity in strongly interacting quan- 
tum field theories from black hole physics," Phys. Rev. Lett. 94, 111601 (2005) 
[arXiv:hep-th/0405231] . 

P. Kovtun and A. Starinets, "Thermal spectral functions of strongly coupled 
N = 4 supersymmetric Yang-Mills theory," Phys. Rev. Lett. 96, 131601 (2006) 
[arXiv:hep-th/0602059] . 

H. Liu, K. Rajagopal and U. A. Wiedemann, "Calculating the jet quenching parameter 
from AdS/CFT," |arXiv:hep-ph/0605178[ 

C. P. Herzog, A. Karch, P. Kovtun, C. Kozcaz and L. G. Yaffe, "Energy loss of a heavy 
quark moving through N = 4 supersymmetric Yang-Mills plasma," JHEP 0607, 013 
(2006) [arXiv:hep-th/0605158| . 

A. Buchel, "On jet quenching parameters in strongly coupled non-conformal gauge the- 
ories," |arXiv:hep-th/0605178| 

S. S. Gubser, "Drag force in AdS/CFT," |arXiv:hep-th/0605182[ 

S. Caron-Huot, P. Kovtun, G. D. Moore, A. Starinets and L. G. Yaffe, "Photon and 



dilepton production in supersymmetric Yang-Mills plasma," arXiv:hep-th/0607237 



P. Binetruy, C. Deffayet, U. Ellwanger and D. Langlois, "Brane cosmological evo- 
lution in a bulk with cosmological constant," Phys. Lett. B 477, 285 (2000) 
[arXiv:hep-th/9910219] . 

E. Kiritsis, "Holography and brane-bulk energy exchange," JCAP 0510, 014 (2005) 
[arXiv:hep-th/0504219] . 

M. Cvetic, S. Nojiri and S. D. Odintsov, Phys. Rev. D 69, 023513 (2004) 
[arXiv:hep-th/0306031] . 



9 



[18] For solving gravity equations we have used the Mathematica code diffgeo5.m by M. 



Headrick, http:/ /www. stanford.edu/~headrick/ 



[19] S. de Haro, S. N. Solodukhin and K. Skenderis, "Holographic reconstruction of spacetime 
and renormalization in the AdS/CFT correspondence," Commun. Math. Phys. 217, 595 
(2001) [arXiv:hep-th/0002230| . 

[20] M. Henningson and K. Skenderis, "The holographic Weyl anomaly," JHEP 9807, 023 
(1998) | arXiv:hep-th/9806087| . 

[21] A. Ashtekar and B. Krishnan, "Dynamical horizons and their properties," Phys. Rev. D 
68, 104030 (2003) garXiv:gr-qc /0308033]. 



10 



